Super-sensitive phase estimation with coherent boosted light using parity measurements
Xu Lan1, 2, 3, Tan Qing-Shou1, †
College of Physics and Electronic Engineering, Hainan Normal University, Haikou 571158, China
School of Mathematics and Computational Sciences, Hunan First Normal University, Changsha 410205, China
Key Laboratory of Low-Dimensional Quantum Structures and Quantum Control of Ministry of Education, and Department of Physics, Hunan Normal University, Changsha 410081, China

 

† Corresponding author. E-mail: qingshoutan@163.com

Project supported by the National Natural Science Foundation of China (Grant No. 11665010), the Key Laboratory of Low-Dimensional Quantum Structures and Quantum Control of Ministry of Education, China (Grant No. QSQC1414), and the Scientific Research Fund of Hunan Provincial Education Department, China (Grant No. 17B055).

Abstract

We consider a passive and active hybrid interferometer for phase estimation, which can reach the sub-shot-noise limit in phase sensitivity with only the cheapest coherent sources. This scheme is formed by adding an optical parametric amplifier before a Mach–Zehnder interferometer. It is shown that our hybrid protocol can obtain a better quantum Cramer–Rao bound than the pure active (e.g., SU(1,1)) interferometer, and this precision can be reached by implementing the parity measurements. Furthermore, we also draw a detailed comparison between our scheme and the scheme suggested by Caves [Phys. Rev. D 23 1693 (1981)], and it is found that the optimal phase sensitivity gain obtained in our scheme is always larger than that in Caves’ scheme.

1. Introduction

Interferometers, as an extremely useful and flexible precise measuring tool, play a key role in the field of quantum metrology.[121] There are mainly two classes of interferometers:[2,3] passive (e.g., Mach–Zehnder interferometer (MZI)) and active (e.g., SU(1,1)) ones. The ultimate goal of both setups is to obtain the sub-shot-noise limit for phase estimation. The sensitivity of these setups crucially depends on the input states as well as the detection scheme.

It is well known that, for the MZI to beat the shot-noise limit (SNL), i.e., with N the total photon number, nonclassicality of the input states, such as the single mode squeezed vacuum states[1,22,23] and the two-mode squeezed vacuum states,[2,24,25] is necessary. In contrast to the MZI, the potential advantage of the SU(1,1) interferometer is that even for bright classical coherent sources of the input states, the SNL can also be surpassed.[4,5] In Refs. [4] and [5], the authors have shown that when inputting coherent states into this interferometer the SNL can be beaten through a photon intensity measurement. However, in general, the photon intensity measurement is not optimal, so the best phase estimation precision given by the quantum Cramer–Rao bound (CRB),[2629] that is with F the quantum Fisher information (QFI), is usually hard to be reached in this case.

Recently, parity detection[24,25,3036] has been considered to be a very useful measurement method in MZI with practically many input states, which can achieve the CRB for phase sensitivity. References [24] and [25] have shown that the extreme sensitivity of parity detection is equivalent to measure all the phase-bearing terms, including the off-diagonal terms, in the density matrix of the two-mode light field before the second beam splitter.[24] However, this detection method cannot be used on the SU(1,1) interferometer directly when the beam splitters in MZI are replaced by the nonlinear optical elements. To the best of our knowledge, the optimal measurement scheme for the SU(1,1) interferometer remains absent.

In this work, we combine the advantages of passive and active interferometers by adding an optical parametric amplifier (OPA) before an MZI. We show that our protocol can obtain a better CRB than the pure active interferometer; and this precision can be achieved by a parity measurement. Comparing with other quantum metrology schemes utilizing the parity detection, the main advantage of our scheme is that it has the potential to reach the phase uncertainties far below the SNL, and only needs the cheapest coherent sources. To demonstrate the advantage of our scheme, we also compare the optimal phase sensitivity gains (OPSG) with respect to the SNL between these two schemes. It is found that the gain obtained in our scheme is always higher than that obtained in Ref. [1] in the whole range of r.

2. Phase estimation in a hybrid interferometer

We consider a passive and active hybrid interferometer, which is formed by an OPA, two 50:50 beam splitters, and two phase shifters, as shown in Fig. 1(a). This setup requires a strong pump beam in addition to two coherent light beams labeled by |α0〉 and |β0〉. The OPA performs the unitary transformation S(r) = exp [r(abab)], where r is the two-mode squeezing parameter, which is relevant to the strength of the nonlinearity and the intensity of the pump beam. The beam splitter is described by ; and the phase-shift unitary operator is U = exp[i (φaaa + φbbb)] = exp(iNdϕd/2)exp(iNsϕs/2), where we have introduced the difference and sum phase shifts ϕd = φaφb and ϕs = φa + φb, which can be precisely estimated by the setups in Figs.1(a) and 1(b), respectively. The corresponding difference and sum number operators are given by Nd = aabb and Ns = aa + bb, respectively.

Fig. 1. (color online) (a) Schematic diagram of the passive and active hybrid interferometry, which is formed by adding an OPA before an MZI. This setup is sensitive to the differential phase shift ϕd = φaφb, and the information is extracted by a parity detection on one of the second beam splitter’s outputs. (b) Schematic drawing of an ideal SU(1,1) interferometer. The beam splitters of the MZI have been replaced by OPAs. There is a π phase shift between the two OPAs. Different from the setup in panel (a), the SU(1,1) interferometer is sensitive to the sum phase shift ϕs = φa + φb and the total photon intensity measurement is taken on the second OPA’s outputs.

Assuming there is no loss in the configuration, then for the setup shown in Fig. 1(a), we can obtain the QFI F = Nc cosh (4r) + 2|α0β0|sinh(4r) + sinh2(2r), where Nc = |α0|2 + |β0|2 is the total amount of coherent light with α0 = |α0|eiθ1 and β0 = |β0|eiθ2. To obtain the optimal QFI, here we consider the phase-matching condition and choose the phases θ1 = θ2 = 0. Therefore, the best phase sensitivity given by the CRB, i.e., , for our scheme is which is far below the SNL, i.e., . Here NT = Nc cosh (2r) + ns is the total photon number inside the interferometer, where the first term corresponds to contributions from the stimulated parametric process, and the second term results from the spontaneous process with ns = 2 sinh2r. In the limit of |α0| = |β0| → 0, equation (1) reduces to the case of two-mode squeezed vacuum states[2,24,25] and the CRB is , which realizes the Heisenberg limit.

Following the same procedure, we can also obtain the phase sensitivity of the pure SU(1,1) interferometer (Fig. 1(b)) for the same input states. When setting θ1 = θ2 = 0, the CRB turns out to be Comparing Eq. (1) with Eq. (2), we can easily find that the quantum CRB obtained in our hybrid interferometer is better than that obtained in the pure SU(1,1) interferometer, since cosh (4r) > sinh2(2r). Moreover, as we know that the CRB can very well quantify the potential advantage of a certain scheme in quantum metrology, but not every detection method is able to achieve this advantage. For example, the precision given in Eq. (2) cannot be realized through the current measurement scheme, such as photon intensity detection and parity detection. In Refs. [4] and [5], the authors obtained the phase sensitivity in some certain cases by detecting the total photon intensity on the second OPA’s outputs of the SU(1,1) interferometry. It is clear to see that this precision does not reach the CRB. Moreover, we can check that for the SU(1,1) interferometry, the parity detection scheme is even worse than the photon intensity detection when |α0|,|β0| > 0.

3. Achieve the quantum CRB by parity detection

Next, we will demonstrate that, in our hybrid scheme (Fig. 1(a)), the sensitivity in Eq. (1) can be obtained by just implementing the parity detection, which was first introduced into optical interferometry by Gerry,[31] and has been applied to precise phase estimation. When specifying a as the mode to measure, the signal of the parity detection is 〈Πa〉 = 〈(−1)aa〉.

In our scheme, the expectation value of the parity operator 〈Πa〉 in mode a can be obtained as where Note that, in Eq. (3) the parity signal 〈Πa〉 does not depend on the sum phase shift ϕs = φa + φb.

The phase sensitivity Δϕd of an interferometer is characterized by using the error propagation formula Δϕd = ΔΠa/|Πa〉/∂ϕ〉, where . When |α0|2 = |β0|2 = Nc/2, the explicit form of the optimal phase sensitivity with parity detection reads in the limit ϕd → −π/2. Notably, equation (4) saturates the quantum CRB given in Eq. (1) since 2|α0β0| sinh (4r) = Nc sinh(4r) in this case. This implies that the parity detection is the optimal measurement scheme for our hybrid interferometer.

In Fig. 2, we plot the parity signal 〈Πa〉 and the phase sensitivity Δϕd as a function of differential phase ϕd = φaφb. There is a central peak at ϕd → −π/2 for the parity signal, and the width of the peak becomes narrower and narrower when increasing the value of r (see Fig. 2(a)). This leads to an increased magnitude of slope (i.e., increased |Πa〉/∂ϕ|) around ϕd → −π/2, and hence the phase sensitivity is improved. From Fig. 2(b), we can clearly see that the best phase sensitivity given by CRB is found at ϕd → −π/2; and for fixed photon number Nc, when increasing r the phase uncertainty decreases gradually. This point can be seen from Eq. (4).

Fig. 2. (color online) (a) The parity signal 〈Πa〉 and (b) the phase sensitivity Δϕd as a function of differential phase ϕd = φaφb when r takes different values, where |α0|2 = |β0|2 = Nc/2 = 10.

We now use the numerical techniques to investigate the optimal phase sensitivity Δϕdmin in the case of |α0| ≠ |β0|, namely, k = |α0|2/Nc, by setting Nc = 500. From Fig. 3, we find that the optimal sensitivity is reached around k = 0.5, but it is also shown that, if r > 1.5 the values of the phase sensitivity hardly have changes when Nc = 500. In other words, it is insensitive to the value of k. Hence, we can safely say that as long as r is not too small, for a fixed total amount of coherent light Nc, the phase sensitivity given in Eq. (4) can still be acquired.

Fig. 3. (color online) The quantum CRB Δϕdmin as a function of the fraction k ≡ |α0|2/Nc when r takes different values. Here the total coherent photon number is Nc = 500.

It is worthwhile noting that in the high power regime (i.e., for large Nc or large r), under our consideration, the accurate number-resolving may be a difficult task. However, there are some proposals[25,30] which can implement the parity detection while they do not need to have photon number counting. These efforts confirm that the high-intensity input states scheme can be constructed with current technology, and then reach the sub-SNL of phase sensitivity. Below, we shall compare our method with another high-intensity quantum-enhanced method presented by Caves.

4. Comparison with Caves’ scheme

The scheme presented by Caves employs a high-intensity coherent state and a low-intensity squeezed vacuum state as the input states of an MZI,[1] and the quantum CRB turns out to be[22,35] . Such a phase sensitivity can be realized either by implementing parity detection[25] or using Bayesian parameter estimation.[22] Now, we will demonstrate the advantages of our scheme when comparing with Caves’ scheme. First, we consider an interesting case of Nc ≫ sinh2r, from Eq. (4) we immediately have , while for Caves’ scheme we have . This means that our scheme is exponentially better than Caves’. However, this is not to say that our scheme is directly more efficient, since the total photon numbers are different between these two schemes. To make a fair comparison, below we introduce the OPSG with respect to the SNL, which is defined as[12] For our scheme, we have Whereas for Caves’ scheme, the OPSG is given by since the total photon number in this scheme is . If , it indicates that our scheme is more efficient than Caves’.

In Fig. 4, we compare the OPSG between our scheme (gmin) and Caves’ scheme ( ) as a function of squeezing parameter r for various coherent photon numbers Nc. As can be seen, the phase sensitivity gains increase when increasing the squeezing parameter r, and we always have in the whole range of r. For our scheme, when r = 1 it is possible to overcome the SNL by 5.8 dB, while reaching 14.6 dB when r = 3 (this two-mode squeezed parameter is available in current experience[37]); and these results are independent of the coherent photon numbers Nc. In contrast to these, for Caves’ scheme we have with r = 1, and for larger r there are bifurcations when increasing coherent photon numbers Nc. It can be seen that the larger Nc is, the better OPSG can be obtained. For example, when r = 3 we have with Nc = 50, while with Nc = 1000.

Fig. 4. (color online) A comparison of the optimal phase sensitivity gains between our scheme gmin and Caves’ scheme as a function of squeezing parameter r for various coherent photon numbers Nc. Note that in our scheme the values of OPSG are independent of Nc.
5. Conclusion

In summary, we have investigated a hybrid scheme for optical interferometry, which assembles the advantages of both passive and active interferometers. This scheme can reach a scaling far below the SNL for phase estimation with bright sources. Comparing with the pure active interferometer schemes with coherent light input, our protocol can obtain a better QCR, and precisely this can be achieved when a parity measurement is utilized. Furthermore, we have compared our scheme with Caves’ scheme; it is found that our scheme is not only more sensitive, but also more efficient when the same total photon number is input.

Reference
[1] CavesC M 1981 Phys. Rev. 23 1693
[2] YurkeB McCallS L KlauderJ R 1986 Phys. Rev. 33 4033
[3] SandersB C MilburnG J ZhangZ 1997 J. Mod. Opt. 44 1309
[4] Plick W N DowlingJ P AgarwalG S 2010 New J. Phys. 12 083014
[5] MarinoA M Corzo TrejoN V LettP D 2012 Phys. Rev. 86 023844
[6] HollandM J BurnettK 1993 Phys. Rev. Lett. 71 1355
[7] DornerU Demkowicz-DobrzanskiR SmithB J LundeenJ S WasilewskiW BanaszekK WalmsleyI A 2009 Phys. Rev. Lett. 102 040403
[8] GiovannettiV LloydS MacconeL 2011 Nat. Photon. 5 222
[9] SandersB C MilburnG J 1995 Phys. Rev. Lett. 75 2944
[10] MaJ WangX SunC P NoriF 2011 Phys. Rep. 509 89
[11] HumphreysP C BarbieriM DattaA WalmsleyI A 2013 Phys. Rev. Lett. 111 070403
[12] PezzeL SmerziA 2013 Phys. Rev. Lett. 110 163604
[13] DowlingJ P 2008 Contemp. Phys. 49 125
[14] LvovskyA I SandersB C TittelW 2009 Nat. Photon. 3 706
[15] EscherB M de Matos FilhoR L DavidovichL 2011 Nat. Phys. 7 406
[16] HudelistF KongJ LiuC JingJ OuZ Y ZhangW P 2014 Nat. Commun. 5 3049
[17] YueJ D ZhangY R FanH 2014 Sci. Rep. 4 5933
[18] TanQ S LiaoJ Q WangX Franco N 2014 Phys. Rev. 89 053822
[19] LiuG Q ZhangY R ChangY C YueJ Do FanH PanX Y 2015 Nat. Commun. 6 6726
[20] WangY T TangJ S HuG WangJ YuS ZhouZ Q ChengZ D XuJ S FangS Z WuQ L LiC F GuoG C 2016 Phys. Rev. Lett. 117 230801
[21] WeiC P HuX Y YuY F ZhangZ M 2016 Chin. Phys. 25 040601
[22] PezzéL SmerziA 2008 Phys. Rev. Lett. 100 073601
[23] JarzynaM Demkowicz-DobrzańskiR 2012 Phys. Rev. 85 011801
[24] AnisimovP M ChiruvelliA PlickW N HuverS D LeeH DowlingJ P 2010 Phys. Rev. Lett. 104 103602
[25] PlickW N AnisimovP M DowlingJ P LeeH AgarwalG S 2010 New J. Phys. 12 113025
[26] CramerH 1946Mathematical Methods of StatisticsPrincetonPrinceton University500
[27] HelstromC W 1976Quantum Detection and Estimation TheoryNew YorkAcademic Press20
[28] HolevoA S 1982Probabilistic and Statistical Aspects of Quantum TheoryAmsterdamNorth-Holland110
[29] BraunsteinS L CavesC M 1994 Phys. Rev. Lett. 72 3439
[30] RoyerA 1977 Phys. Rev. 15 449
[31] GerryC C 2000 Phys. Rev. 61 043811
[32] GerryC C MimihJ 2010 Contemp. Phys. 51 497
[33] GerryC C MimihJ 2010 Phys. Rev. 82 013831
[34] ChiruvelliA LeeH 2011 J. Mod. Opt. 58 945
[35] SeshadreesanK P AnisimovP M LeeH DowlingJ P 2011 New J. Phys. 13 083026
[36] SeshadreesanK P KimS DowlingJ P LeeH 2013 Phys. Rev. 87 043833
[37] MartiniF D SciarrinoF VitelliC 2008 Phys. Rev. Lett. 100 253601